What drives the translocation of stiff chains? 
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We study the dynamics of the passage of a stiff chain 
through a pore into a cell containing particles that bind re- 
versibly to it. Using Brownian Molecular Dynamics simula- 
tions we investigate the mean-first-passage time as a function 
of the length of the chain inside, for different concentrations of 
binding particles. As a consequence of the interactions with 
these particles, the chain experiences a net force along its 
length whose calculated value from the simulations accounts 
for the velocity at which it enters the cell. This force can 
in turn be obtained from the solution of a generalized diffu- 
sion equation incorporating an effective Langmuir adsorption 
free energy for the chain plus binding particles. These re- 
sults suggest a role of binding particles in the translocation 
process which is in general quite different from that of a Brow- 
nian ratchet. Furthermore, non-equilibrium effects contribute 
significantly to the dynamics, e.g., the chain often enters the 
cell faster than particle binding can be saturated, resulting in 
a force several times smaller than the equilibrium value. 

I. INTRODUCTION 

The transfer of genetic material through the membrane 
surrounding a cell nucleus is fundamental to the under- 
standing of basic cell processes from gene therapy to viral 
infection. The motion of linear polymers through pores 
into confined volumes also arises in many other biologi- 
cal contexts [1], perhaps the most common examples of 
which include the translocation of proteins from the cy- 
tosol into the endoplasmic reticulum, or into mitochon- 
dria or chloroplasts. The export of m-RNA through nu- 
clear pore complexes is still another example of great 
importance. Despite the longstanding and widespread 
interest in this process, however, our knowledge about it 
is still rudimentary. 

The process of translocation under influence of an ex- 
ternal field or a chemical potential gradient has recently 
been studied extensively [2-5] . There have also been sev- 
eral theoretical studies which have specifically investi- 
gated chain translocation in the presence of binding par- 
ticles [6-9]. The major role of these binding particles 
has been recognized as a "Brownian ratchet", a mech- 
anism which was introduced over ten years ago in pi- 
oneering work by Simon, Pcskin and Oster [6]. They 
happen to have treated the case of protein translocation, 
but their arguments apply equally well to nucleic acids. 
To account for translocation rates fast compared to sim- 
ple diffusion under a wide variety of conditions and cir- 
cumstances, they proposed that non-specific binding by 



globular proteins results in a "biased" — or "ratcheted" 
— motion of the chain. Many experiments confirm that 
efficient translocation can indeed take place without the 
involvement of motor proteins. For example, the entire 
length (about 40 microns) of the DNA which comprises 
the genome of T5 phage is observed to enter its bacterial 
cell host without requiring metabolic energy [10]. The 
experiment of Salman et al. on phage A DNA shows sim- 
ilarly that translocation of a comparable length of chain 
can occur without the help of active processes; simple 
diffusion would require significantly longer times [11]. 

Each of the mechanisms mentioned above for chain 
translocation, namely diffusion and ratcheting, corre- 
sponds to a different time scale and to different physics. 
Simple diffusion requires a characteristic time td — 
L 2 /2D, where L is the total length of the polymer and 
D its diffusion coefficient. In the "ratcheting" scenario, 
as soon as a specific length — 6 — of the polymer enters, 
a protein binds to it and the chain is no longer able to 
diffuse backwards because the pore size is too small for 
the DNA/protein complex to pass through. In this case, 
the chain simply diffuses from one binding site to the 
next, and the translocation time is equal to the product 
of the time it takes for the chain to diffuse the distance 
S times the number of ratcheting sites M = L/S, i.e., 
Uatchet = MS 2 /2D = L 2 /(2MD) = t d /M, correspond- 
ing to a speed-up of the translocation time by a factor of 
M over simple diffusion. As pointed out by Simon, Pe- 
skin and Oster, this time represents an idealized limit in 
which the ratcheting mechanism functions " perfectly" , 
i.e., as each successive binding site passes into the cell 
it is bound irreversibly by a protein that prohibits the 
chain from diffusing backwards. In actuality, however, 
an entering site is not necessarily bound immediately by 
a protein, and/or the protein does not stay adsorbed long 
enough to act as a ratchet at that site; accordingly, the 
translocation time is increased beyond LS/2D by a fac- 
tor that depends on the ratio of on- and off- rates for 
binding. 

In this paper we consider explicitly the effect of bind- 
ing particles and show that — via a new mechanism — 
translocation can occur at rates significantly faster than 
that provided by the "perfect" ratcheting scenario de- 
scribed above. More specifically, we argue that the par- 
ticles which bind reversibly to the chain give rise to a net 
force on the chain that pulls it into the cell. Further- 
more, this force accounts fully for the translocation pro- 
cess and embraces the different mechanisms mentioned 
above, e.g., in a special limit the Brownian ratcheting 
appears as a particular idealization (and appealing sim- 
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plification) of the effect of such a force. The magnitude of 
the force depends in a delicate way on the concentration 
of the binding particles and on their diffusion coefficient 
relative to that of the chain. In the overdampcd limit 
the translocation time in the presence of a force F is ex- 
pected to take the form tp = L/v = L(/F, where £ is 
the friction coefficient of the chain, related to its diffusion 
constant through the Einstein relation D = ksT/C,. 

Using Brownian Molecular Dynamics (BMD) simu- 
lations, we calculate separately the average force and 
the translocation time in the presence of binding par- 
ticles and find that they do indeed obey the relation 
tp = L(/F. Consequently, the translocation process is 
force-driven and the translocation time turns out to be 
longer or shorter than the ideal ratcheting time, depend- 
ing on the concentration and diffusion coefficient of the 
binding particles. These results can be understood in 
the context of a generalized Fokker-Planck equation for 
the probability that at time t the chain will have entered 
the cell to a distance x and have n particles bound to it. 
The drift terms corresponding to chain entry and parti- 
cle binding are obtained from derivatives with respect to 
entry distance and binding number of a Langmuir adsorp- 
tion free energy for the overall system. 1 We find that the 
binding process involves important non-equilibrium ef- 
fects in general, on which depends the actual value of the 
force pulling the chain in, and hence the ratio tp/t ratc het- 

II. SIMULATION 

Brownian Molecular Dynamics simulations of a stiff 
polymer translocating through a pore and into a cell 
are performed using a coarse-grained model in which the 
chain is represented by a rigid rod of beads. The beads 
arc rigidly linked to their nearest neighbors along the 
chain and do not interact with each other. In this way 
we model DNA as a perfectly straight — rather than 
the usual semi-flexible — chain, because the focus of our 
work is on the entering segment of chain which is within 
a persistence length from the pore. The link between ad- 
jacent beads is rigid in order to avoid chain contraction 
and extension resulting from the binding of particles; we 
also neglect changes in the shape or twist of the chain 
due to binding. The distance between monomers along 
the chain, a, corresponds to the "footprint" of binding 



particles in that the center of each bead is considered as 
an absorbing site. The binding particles are modeled as 
spherical, interacting with each other through the repul- 
sive part of a Lennard Jones (LJ) potential with diameter 
a; the interaction between binding particles and chain 
monomers is treated by a full (12-6) LJ potential. Since 
the distance between the absorbing sites is equal to the 
diameter of the binding particles, each site can be sur- 
rounded by a maximum of six particles. The cell-particle 
interactions are taken to vanish for particles within the 
radius R s of cell and to increase as (R — R s ) 4 for parti- 
cles at distances greater than R s from the sphere center. 
This potential is simply a convenient way to describe an 
interior wall. Finally, we treat the pore itself as being 
completely " inert" , having no effect on the chain except 
to allow it to enter or leave the cell. 

We focus on the dynamics of the chain once one end 
has been inserted. Let x denote the length of chain inside 
the cell and r t the position of the i th binding particle 
(see Fig. 1). The time evolution of these coordinates is 
described by the overdamped Langevin equations 2 

dr 

= fiD /(k B T) + hi (1) 

^ = FD rod /(k B T) + B . (2) 

Here fi and F are the deterministic forces acting on each 
particle and the rod, respectively, and bi and B are the 
corresponding random — "Brownian" — forces satisfy- 
ing < bi(i) >= 0, < bi(i) • bj(t') >= 6D 5{t - t')<%, 
< B(t) >= and < B(t)B(t') >= 2D rod S(t-t') >. D is 
the diffusion coefficient of an individual binding particle, 
related to its friction coefficient, ( , through the Einstein 
relation D — k B T/Q . As for the chain, we introduce 
an effective diffusion coefficient, D roc [ — kBT/C, rod which 
in principle may include all the pore-DNA interactions. 
Since so little is known about these complicated interac- 
tions, we have simply taken D rod = Dq / (L / a) , consistent 
with the translational diffusion coefficient of a stiff chain 
being inversely proportional to its length. In all that fol- 
lows we use a, ksT, and cr 2 /D as the units of length, 
energy, and time, e, the LJ binding energy between the 
particles and the monomers, is set equal to 5fc^T; the 
diameter of the spherical cell is 2R S = 24cr and the total 
length of the chain is L = 16a. 



1 In the simplest case of the fully equilibrated, fixed chemical 
potential Langmuir adsorption problem, the ID pressure of 
the system corresponds to the force on the chain: Pm = f = 
(fcsT/5)ln(l + 0exp(/?e)) where <f> is the volume fraction of 
binding particles, e their binding energy, and 8 the binding 
site size. In the limit of large binding energies (/3e S> 1) one 
recovers the / ~ e/8 result discussed in Section IV; otherwise, 
entropic effects are important as well. We thank Pierre-Gilles 
de Gennes for these observations. 



2 We ignore the force in the y and z directions on the chain, 
and let the rod move only in the x direction. 
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FIG. 1. A schematic view of the translocation of the rigid 
chain into the spherical cell. 

Figure 1 is a schematic snapshot of the simulated sys- 
tem. The particles bind predominantly at the tip and 
then move back along the chain, with particles occasion- 
ally adding in empty spaces along the rod; at the same 
time, particles leave from other parts of the chain and 
allow for new particles to bind. Geometrically, up to six 
particles can sit around each monomer, but just as in the 
familiar Langmuir adsorption problem there are always 
empty sites on the chain due to entropic factors. The fact 
that particles are mostly added at the tip is completely 
a dynamical effect. Under the influence of particle bind- 
ing, the chain is in general moving too fast to allow for 
saturation of adsorption along its length. This results 
in less than the equilibrium number of occupied sites on 
the chain. In general, particles bind to the tip of the 
chain by pushing aside some already attached particles. 
However, the force pulling in the chain is exerted mostly 
at the entering positions where empty binding sites first 
appear. 

III. RESULTS 

Figure 2 shows the Mean First Passage Time (MFPT) 
versus the length of the chain inside the cell, x. MFPT 
is the average, over a large number of trajectories, of the 
time it takes for the front tip of the chain to first arrive 
to the position x. Each curve in the figure corresponds 
to a different value of N, the number of binding particles 
inside the cell. 

As mentioned in the INTRODUCTION, there exist 
different mechanisms for the translocation of the rod un- 
der the influence of the binding particles. If the chain 
simply diffused into the spherical cell, the MFPT would 
be equal to td — L 2 /(2D ro d) = 2048er 2 /D correspond- 
ing to the length L — 16a; see the quadratic function 
depicted by the heavy solid curve in Fig. 2. As shown 
in the figure, td lies significantly above the translocation 
times we find in our simulations; even in the presence 
of only five binding particles, the MFPT is about three 
times shorter. 
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FIG. 2. Mean First Passage Time (MFPT) as a function 
of entry distance x, for each of several different values of the 
number of binding particles N, as calculated from our BMD 
simulations. The thick solid curve describes the MFPT vs 
x for simple diffusion of the chain into the cell. The dotted 
curve shows the MFPT for the case where translocation would 
occur via perfect ratcheting (see text). 

The dotted line in Fig. 2 represents the time that it 
would take for the entire chain to enter if there were 
ratchets functioning perfectly at every absorbing site. In 
this case, the chain simply undergoes successive and in- 
dependent diffusions between neighboring sites, complet- 
ing each in a time a 2 / (2D ro d); the MFPT is equal to 
the product of this time and the number of steps (x/a) 
associated with the entry distance x. According to the 
ratcheting mechanism, then, the slope of time versus dis- 
placement is simply a/(2D ro d) =8a/Do, for D ro d = 

The numerical results from our simulations, as well 
as the theory outlined in the next section, confirm the 
presence of a quite different translocation mechanism, 
namely, drift due to a net force exerted on the chain 
by binding particles. Entry into the cell of successively 
longer portions of chain "feeds" new binding sites into 
the system; as each additional particle binds to the chain 
the free energy of the chain drops, and this reduction 
gives rise to a force pulling the chain into the cell. In 
the presence of a constant force, and in the overdamped 
limit, the corresponding translocation time is tp = L/v 
where v = FD ro d/kBT is the velocity of the chain. 

With competing mechanisms operative, the MFPT will 
reflect predominantly the one with the smallest translo- 
cation time. As already remarked, the diffusion time is 
always much longer than those arising from the other two 
mechanisms (see Fig. 2). However, the ratcheting time 
could be longer or shorter than that of the force-driven 
process according to whether the deterministic force F is 
larger or smaller than the "effective" Brownian ratcheting 
force. Comparison of the ratcheting velocity of a chain, 
2D ro d/<7, with the usual expression for the velocity of 
a chain under a constant force, v = F ratc hetD ro d/ ksT , 
shows that the "effective" Brownian ratcheting force is 
2ksT / a. Consequently, the ratcheting mechanism is ex- 
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pected to be dominant in the presence of weak enough 
driving forces. Nevertheless, even in this limit, as we shall 
see below, the attractive force is operative and completely 
determines the translocation velocity. 

Several important observations can be extracted from 
Fig. 2. The first is that translocation times depend on 
the number (concentration) of binding particles. This is 
also true, of course, for the ratcheting mechanism, since 
higher concentrations of binding particles imply faster 
"on-rates" (fc+), hence higher translocation velocities. 
But, as Simon, Peskin and Oster have themselves em- 
phasized, there is a maximum translocation rate corre- 
sponding to the limit of large This maximum 
rate — or minimum time, t ratc het = L5/2D ro d — also 
corresponds to the smallest distance (a in the present 
model) between ratcheting (binding) sites. From Fig. 2 
we see that the translocation occurs even faster than the 
limiting ratcheting prediction when the number of bind- 
ing particles exceeds N = 100. This result suggests that 
an additional mechanism is operative, which we show is 
associated with a net force acting on the chain along its 
direction of motion. 

More explicitly, the slope of MFPT vs x plots reveals 
that the average velocity of the chain, v, remains almost 
constant throughout the translocation process (except 
right at the beginning and towards the end). From the 
relation F = CrodV one expects that the average force 
on the chain also stays constant during this process, in 
which case the slope of time t versus < x > will be almost 
the same as that of MFPT versus x iff vx/D ro d ^ l 3 . 
Using F = Crod,v, we calculated the effective force on the 
chain from the velocity obtained in our MFPT vs length 
curves (Fig. 2). Specifically, we calculated the average 
velocity of the chain over the range x = 4cr to x = llcr, 
i.e., in the region where the MFPT versus x curves are 
almost linear (see Fig. 2). 

Alternatively, we can also determine the force acting on 
the chain directly from the simulations. The squares in 
Fig. 3 represent the average force exerted on the chain by 
the binding particles (coarse-grained over er), as a func- 
tion of the length inside. This force is seen to be nearly 
constant; furthermore, it is found to agree within nu- 
merical uncertainty, as shown in Fig. 4, with the force 
calculated (as described above) from the average velocity 
(i.e., the inverse slope in Fig. 2) for all the values of N 
(up to 300) that we treated. 
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FIG. 3. The squares illustrate the force calculated directly 
in the simulation (coarse-grained overer), as a function of the 
length of the chain inside, for the case N = 100. The circles 
show the force calculated from an identical BMD simulation 
but with a rod diffusion coefficient which is 16 times smaller. 
The dashed curve corresponds to the force derived from solu- 
tion of the coupled entry /adsorption diffusion equation in the 
approximation of quasi-equilibrium. The solid curve is the 
solution of the full Fokker-Planck Eq. (4) using a = 10.8<r 
(see text). 

For N < 50, where the forces are of the order of unity 
or smaller (i.e., smaller than the "effective" Brownian 
ratchet force 2ksT/a), Fig. 2 shows that the translo- 
cation times begin to be significantly smaller than the 
perfect ratcheting limit. This scenario is in principle 
embodied in the basic result of Ref. [6] in which the 
translocation time is written as t ra tchet(l + 2K~ 1 ) where 
tratchet — L5/2D ro d is the ideal ratcheting time and 
K = k + /k- is the effective strength of binding, expressed 
as a ratio of "on" and "off" rates. K 1 corresponds 
to saturated binding and to ideal ratcheting; otherwise, 
the ratchet mechanism becomes less efficient and the 
translocation times are longer than t ratc het ■ ft is notable 
that even in this weak force situation, where one expects 
the ratcheting mechanism to dominate, we still observe 
a force-controlled translocation (i.e., the translocation 
time is still determined by the force pulling the chain in, 
as shown in Fig. 4). The reason is that these weaker 
forces correspond to smaller numbers of less strongly 
bound particles, and hence also to less efficient ratch- 
eting. 



3 The translocation time of a polymer can be defined in 
terms of the mean first passage time r(x), which satisfies the 
backward Fokker-Planck equation with a reflecting bound- 
ary condition at the hole and an absorbing boundary con- 
dition at x. The solution to the backward Fokker-Planck 
equation in the presence of a constant force F = fi> is 
t(x) = B + £(e—/ c - 1). If v/D » 1, then r(x) = §. In 

the other limiting case, i.e., when v/D <C 1, then t(x) = 
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FIG. 4. Comparison of the force obtained directly in the 
simulation (hollow circles) with the one calculated using the 
relation F = ( ro dV (filled circles) where v is the velocity cor- 
responding to the inverse slope of MFPT vs. x data. 



IV. DYNAMICAL THEORY 

The translocation process simulated above can be de- 
scribed theoretically using a simplified dynamical model. 
The two relevant variables are the length of the chain x 
inside the cell and the number of particles n attached to 
the chain. As the chain enters a length x, the number of 
available absorbing sites increases. With e the binding 
energy of a single particle we can write the (Langmuir 
adsorption) free energy, A, of the system as 

(6x/cr)! 



A(x, n) — —ne — kgT log 



n\(6x/a — n)\ 
k B T{N - n) log ■ ' 



(3) 



(N-n)v ' 

Here N is the total number of particles, V is the volume 
of the spherical cell, and Vq is the volume of a single 
particle. The coefficient 6 in the second term of Eq. (3) 
is the number of particles that can interact attractively 
with a single binding site (chain monomer); accordingly, 
6x/cr is the total number of available sites on a chain of 
length x. The first term in Eq. (3) is the energy gain 
due to binding; the second is the entropic contribution 
associated with the total number of ways in which 6x/a 
sites can be occupied by n particles; and the last term 
is the (ideal gas) contribution associated with the "free" 
particles, numbering N — n. 

We consider the dynamics of translocation as a coupled 
diffusion process involving both the x and n degrees of 
freedom. Using Mesoscopic Nonequilibrium Thermody- 
namics (MNET) [12], one can derive the Fokker-Planck 
equation governing the time-dependent probability den- 
sity p(x, n, t) that at time t a length x of the chain has 
passed through the hole and has n particles attached to 
it: 



dp(x, n, t) 
di 



d_ 

dx 
d 



dn n \ ksT 



1 dA(x,n) 
ksT dx 
( 1 dA(x, n) 



P 



dn 



dp 
dx 

dp 

dn 



(4) 



Note that there is a drift and a diffusion term for each of 
the x and n variables. The force driving the translocation 
is F(x,n) = — dA ^ x ' n ' > , and we see from Eq. (3) that its 
origin is entropic, arising from the second ("Langmuir") 
term in Eq. (3); the binding of the particles gives rise 
to a force that pulls the chain in. Similarly, the factor 
— dA g^ n ^ is a "thermodynamic force" driving the particle 
binding. D ro d is the spatial diffusion coefficient of the 
rod (with the usual dimensions of length 2 time~ 1 ), while 
D n (which in general may depend on x and n) is the 
kinetic rate constant (with dimensions of time -1 ) for the 
process of particle binding and unbinding. A crude, but 
time-honored and physically reasonable, expression for 
D n comes from the Smoluchowski theory of aggregation 
dynamics [13] for diffusing particles: 



N 

D n = a— — Ah 



V, 



cell 



(•5) 



where a is a length of order the particle size. D n here 
is simply proportional to the concentration of binding 
particles N/V, and their spatial diffusion coefficient Do- 
The Fokker-Planck Eq. (4) provides a complete de- 
scription of the kinetics of both chain entry and particle 
binding. However, a simpler description can be achieved 
by considering the possibility of time scale separation. 
The characteristic times for the entry and the binding 
processes scale as t x ~ l/D roc i and t„ ~ 1/D n , respec- 
tively. If the binding process is very fast compared to the 
chain entry (i.e. T n / t x ~ D ro( i/D n <C 1), it is reasonable 
to assume that the "fast" variable — here the number 
of attached particles, n — will decay very rapidly to its 
equilibrium distribution. In this case the process can be 
described by the evolution of the slow variable, the po- 
sition x of the chain. Suppression of the fast variable 
can be carried out using the standard adiabatic elimi- 
nation technique [14], which is essentially equivalent to 
integration of the Fokker-Planck equation over the equi- 
librium distribution of the fast variable. The resulting 
one dimensional Fokker-Planck equation is 



dp(x,t) = d_ D 
dt dx r 



F(x) dp 
ksT dx 



where 



. dA(x,n) 
^ / ^ f eq (n;x)dn, 



(6) 



(7) 



is the average driving force. Here we have assumed that 
the spatial diffusion coefficient of the rod is indepen- 
dent of n, and defined p(n,x,t) = p(x,t)f eq (n;x), with 
f eq (n;x) the "local" equilibrium distribution, 



f eq (n;x) — AQ.n) 

J e k B T dn 



(8) 



In the particular case where f eq (n; x) = 5(n — n eq (x)), 
the force driving the translocation process becomes 
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F eq (x) = 



_ dAjx.n) 



dx 



and the number of attached par- 



ticles is equal to the equilibrium one, n eq (x), given by the 



solution of dA Q^ 



0. 



The above quasi- (or "local") equilibrium approach 
basically assumes that, as soon as the chain's advance 
makes available new sites, particles bind to them im- 
mediately. Figure 3 compares, for N = 100, the av- 
erage force calculated in our simulation (squares) with 
the "quasi-equilibrium" force given by Eq. (7) (dashed 
curve). The comparison shows that the actual force that 
the binding particles exert on the chain is significantly 
smaller than the one which follows from the assumption 
that binding equilibration can keep up with chain entry 
(and this is true for all other N values considered). To 
check the limit in which the simple quasi-equilibrium ID 
Fokker-Planck description, (6), provides an accurate de- 
scription of the translocation process, we slowed down 
the entry of the rod by decreasing D ro d- The circles in 
Fig. 3 shows the average force calculated via BMD sim- 
ulation when D ro d = 1/960 (vs 1/16). Agreement with 
the "quasi-equilibrium" force evaluated from (7) is now 
excellent, confirming the validity of our model for the 
thermodynamic free energy A(x,n). Comparably good 
agreement is found between the number of attached par- 
ticles calculated in the simulation and that predicted by 
the quasi-equilibrium theory. (Alternatively, we can in- 
crease the diffusion coefficient D of the binding parti- 
cles, in which case the calculated force also approaches 
its quasi-equilibrium value, as we checked in the simula- 
tion). The free energy given in Eq. (3) was written in the 
continuum limit under the assumption that the system 
is dilute. Considering its simplicity, the accuracy of the 
model in describing the system is surprisingly good. 

In general, to compare simulation results with the pre- 
dictions of our theory, we need to solve the full Fokker- 
Planck equation for the coupled x and n degrees of free- 
dom. The numerical solution is obtained by converting 
Eq. (4) into its equivalent set of Langevin equations [14], 
and then solving these equations using standard stochas- 
tic algorithms. From these numerical solutions the aver- 
age force, the mean first passage time, and the average 
number of adsorbed monomers can be determined. Us- 
ing the approximate expression for D n , given by Eq. (5), 
with a = 10.8cr (which in fact is very close to the value 
a = Alter from the Smoluchowski theory) we have solved 
Eq. (4) for different values of N. The resulting MFPT 
values vs x, and averaged pulling forces, compare well 
with those from the simulations. As an example, the re- 
sulting average force for N = 100 is represented by the 
solid line in Fig. 3, showing an excellent agreement with 
the results of the simulation. 



V. CONCLUSION 

The present work has attempted to provide a basic 
theoretical framework for treating the translocation dy- 
namics of a stiff chain as it moves into a cell containing 
particles which interact attractively with it and bind to 
it. We conclude that this process is in general force- 
controlled. Obviously, pure diffusion occurs only in the 
absence of any binding particles; and the rectification of 
diffusion — which is the essence of the Brownian ratchet- 
ing mechanism — appears as a manifestation of the bind- 
ing force in a special limit of interaction potentials and 
particle concentration. 

To examine the extent to which the simple Brownian 
ratchet mechanism can account for chain translocation 
due to binding particles, we have performed several sim- 
ulations using different values for S, the distance between 
binding sites. We varied 5 from la to 4a, i.e., from S 
smaller than to larger than the range of the Lennard- 
Jones attractive interaction (about 2a) . When S is suf- 
ficiently larger than the range of interaction — say for 
5 = 4a — we find that the chain performs a cycle of drift 
and diffusion movements. As an absorbing site enters 
the cell the chain is first pulled by the attractive force 
acting on the binding site, and then diffuses between ad- 
sorbing sites, during which time no net force is acting on 
the chain. In this case, the effect of the binding force 
is twofold: it pulls the chain in the region where it is 
acting, and also impedes its backwards diffusion. In the 
limit where the range of the force is small compared to 
the distance between sites, the contribution of the drift to 
the translocation is negligible, but nevertheless the force 
is still rectifying the diffusion. In these circumstances, 
the effect of the force can be most simply described by 
the Brownian ratchet idea. Specifically, the rectification 
arises from the free energy penalty for moving an absorb- 
ing site out of the cell (as provided approximately by our 
Langmuir adsorption model). The efficiency of the rec- 
tification depends on the ratio of free energy penalty to 
the thermal energy fc^T. 

We have also found that the effective force can be 
significantly (several times) smaller than the value one 
would estimate from a quasi-equilibrium treatment of 
the binding and entry dynamics, i.e., from assuming that 
the time scale for particles to diffuse and bind is much 
shorter than that for chain entry. More explicitly, only 
if we increase sufficiently the friction coefficient of the 
chain do we find agreement between the force calcu- 
lated from the simulations and that obtained from the 
quasi-equilibrium solution to the coupled entry/binding 
Fokker-Planck equation. This result indicates that there 
can be significant differences between the stalling force 
measured in single-molecule experiments and the actual 
value of the force during the process of translocation. 

In summary, as established directly from our BMD 
simulations, we find that chain translocation can be un- 
derstood in terms of the pulling force arising from particle 
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binding. Moreover, this mechanism leads to transloca- 
tion times that can become distinctly smaller than ratch- 
eting. By decreasing the diffusion coefficient of the chain 
relative to that of the binding particles, the simulations 
give a still larger effective force pulling the chain along its 
length into the cell. The maximum value of this force, at- 
tainable in the limit N 3> 1 and x — > oo is precisely 6e/<r, 
i.e. the drop in free energy per unit length when all bind- 
ing sites are occupied. These results are nicely confirmed 
by solving directly a coupled Fokker-Planck equation for 
the chain entry and particle binding dynamics. We con- 
clude that physically realistic situations are in general 
more complicated than a ratcheting mechanism in which 
it is assumed that the only effect of particle binding at 
a chain site is to prohibit its diffusing back through the 
pore. 
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